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ABSTRACT 

Transit timing variations (TTVs) are useful to constrain the existence of perturbing 
planets, especially in resonant systems where the variations are strongly enhanced. 
Here we focus on Laplace-resonant three-planet systems, and assume the inner planet 
transits the star. A dynamical study is performed for different masses of the three bod- 
ies, with a special attention to terrestrial planets. We consider a maximal time-span 
of ~ 100 years and discuss the shape of the inner planet TTVs curve. Using frequency 
analysis, we highlight the three periods related to the evolution of the system: two 
periods associated with the Laplace-resonant angle and the third one with the pre- 
cession of the pericenters. These three periods are clearly detected in the TTVs of 
an inner giant planet perturbed by two terrestrial companions. Only two periods are 
detected for a Jupiter-Jupiter-Earth configuration (the ones associated with the giant 
interactions) or for three terrestrial planets (the Laplace periods). However, the latter 
system can be constrained from the inner planet TTVs. We finally remark that the 
TTVs of resonant three or two Jupiter systems mix up, when the period of the Laplace 
resonant angle matches the pericenter precession of the two-body configuration. This 
study highlights the importance of TTVs long-term observational programs for the 
detection of multiple-planet resonant systems. 

Key words: planets and satellites: detection, planetary systems, celestial mechan- 
ics, planets and satellites: dynamical evolution and stability, planets and satellites: 
individual: Kepler-9, Kepler-18, Kepler-30, Gliese 876 



1 INTRODUCTION 

More than 100 multi-planetary systems have been detectecQ- 
At present, at least three extrasolar systems are believed 
to be in a multiply-resonant configuration, similar to the 
Lapl ace resonance of the Jupi ter's Galilean mo ons: HR 8799 
fe.p. iReidemeister et al.ll2009l ) and Gliese 876 (|Rivera et al.l 
I2OI0I ) in a Laplace 4:2:1 configuration, and the KOL500 five- 
body exosystem tha t presents two three-body resonances 
(|Lissauer et al.ll201ll '). 

In formation theories, it has been shown that, at least 
in the case of low-mass planets such as the ones in our Solar 
System, gas-driven migration can f orce the planets to en - 
ter into a multiple-plane t resonance (t Morbidelfi et al.ll2007l ). 
Recently, iLibert fc Tsigani s ( 2011a„ b|) have shown that a 
Laplace three-planet resonance can be achieved for Jupiter- 
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mass planets and, provided that enough eccentricity damp- 
ing is exerted on the planets, such a three-planet system 
could be stable after the dissipation of the gas disc. 

Multiple-planet resonances can act as a phase- 
protection mechanism, especially for systems consisting of 
a lot of massive bodies in a particul arly compact config- 
uration. In their study of Gliese 876, I Rivera et all (I2OIOI ) 
have shown that the Laplace resonance is essential for the 
long-term stability of the system. A lot of KOIs also show 
evidenc e of their proximity to multiple-planet commensura- 
bilities (|Lissauer et al.|[201ll L but it needs confirmation since 
the orbital parameters of the planets are still unknown. 

The TTVs method is a powerful technique to in- 
fer the existence of additional non-transiting planets 
(lMiralda-Escudell2002l : lAgol et all [20051 : iHolman fc Murray! 
l2005h . Indeed a companion (not necessarily transiting) 
planet interacts gravitationally with the transiting planet 
and generates deviations of the transit times with respect 
to the strictly periodic Keplerian case. This is especially the 
case for resonant systems where the variations are strongly 
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Table 1. Frequency analysis of the TTVs curve and the time 
variations of the resonant angles {9i), Laplace angle (0) and argu- 
ments of the pericenter (lj). Top : the resonant two-planet Gliese 
876 system, see Fig.[Tl Bottom : the Laplace-resonant three-planet 
system, see Fig. (2] The main periods are listed by decreasing am- 
plitude of the trigonometric term (ci is the highest amplitude). 
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enhanced. The TTVs t echnique has been applied using 
groun d-based data (e.g. StefFen fc Ago ]ll2005l: Collier et al.l 
I2OO7I]. MOST data fe.g. iMiller-RiccTet al. 200i), HST data 

g. I A — 1 C4-„fr — iH 



(e.g. lAgol fc S tcffcn 20p' 7|) and recently the sp ace missions 



CoR oT (e.g. Bean|,2009irrCsizmadia et al.lboiol ') and Kepler 
(e.g. lBoruckill2010l . l201ll ) . However, deriving the orbital el- 
ements and the mass of t he perturber from TTVs is a diffi- 
cult inverse problem (e.g. lNesvornv fc Morbidellill20od ). and 
different configurations can generate similar TTVs with the 
same dominant period (e.g. lFord fc Howellll2007 : Boue et "al] 
|2012| ). Despite these difficulties. iNesvornv et alj I (|2012l ) have 
recently achieved the first detection and characterization 
of a non-transiting planet by TTVs for the KOI-872 sys- 
tem. The TTVs of other Kepler confirm ed /candidate sys- 
tems are currently under investigation dFord et al ] l201ll . 
l2012l : ISteffen et all I2OI2I : iFabrycky et"al] I2OI2I ). The TTVs 
method i s also promising for detecting Trojan planets (e.g. 
Ford fc Howel l 20o3) and moons of transiting planets (e.g. 



Kipping||2009l l. 



Here we propose an initial qualitative exploration of the 
theoretical TTVs of a three-planet multiply-resonant sys- 
tem, assuming that only the inner body transits the star. 
We analyze the shape of its TTVs curve on a maximal time- 
span of ~ 100 years, in order to identify relevant features 
of the resonant evolution of the whole system. This allows 
us to define a simple criterion to distinguish three-body res- 
onant systems from two-body ones. We only consider here 
the Laplace 4:2:1 resonance, varying the mass values of the 
three planets, but our results could be transposed to another 
three-planet resonance. 

The paper is organized as follows. In Section [2j we ana- 
lyze the TTVs signal of a Laplace-resonant system composed 
of three giant planets, identifying the periods with the time 
variation of the orbital elements. Section [3] focuses on the 
possible detection of terrestrial planets in multiply-resonant 
systems. A discussion is given in Section U while our results 
are summarized in Section (5] 



<X> 



> 
H 
H 




10000 



300 
200 
100 




300000 



200000 



100000 



-100000 



-200000 



-300000 



5000 10000 
TIME (days) 




50 



150 200 250 300 
transit number 



450 



Figure 1. Dynamical evolution of the two-planet resonant 
Gliese 876 system on 40 year s. The parameters o f planets b and 
c are issued from Table 2 of iRivera et al. ] 1 I2OIOI) . Masses of the 
planets are 0.72 and 2.27 Mj, and the orbital periods are 30.1 and 
61.1 days, respectively. From top to bottom: eccentricities, reso- 
nant angles (0; = Ai — 2A2 + LUi, i = 1,2), pericenter arguments, 
and TTVs of the inner planet. 



2 TTVS SIGNAL OF A LAPLACE-RESONANT 
SYSTEM 

We assume in this work that the inner planet of the system 
transits the star, such that our analysis relies on the TTVs 
curve of this planet only. We first recall the features of the 
TTVs signal of a two-planet resonant system. 

An example of a two-planet resonant dynamics is shown 
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Figure 2. Dynamical evolution (on 40 years) of a Laplace-resonant three-planet system consisting of three giant plane t s of 1.5 Mj with 
orbital periods of 20.1, 40.5 and 81.2 days respectively. The parameters are issued from the work of lLibert &: Tsiganij l|2011lJ l. Time is 
given in days. 
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Figure 3. TTVs of the inner planet of the Laplace-resonant 
three-planet system of Figure [2] The time interval is 40 years. 



in Figure [T] where the 2:1 mean-motion resonant evolution 
of the two- planet Gliese 876 sys tem is given for 40 years (pa- 
rameters of lRivera et al.l (|2010l )). The orbital periods of the 
two planets are 30.1 and 6 1.1 d ays, respect i vely. A s stated by 
iHolman fc Murray! (120051 ) and lAgol etall |2TOi), the TTVs 
curve of the inner planet in Figure [1] (bottom panel) shows 



two periods: a short one associated with the resonant angles 
{9i = Ai — 2X2 + LOi, i — 1,2, X being the mean longitude 
and LU the argument of the pericenter), while the longer one 
corresponds to the precession rate of the arguments of the 
pericenters. Their values are given in Table [1] (t op) where, 
by resorting to frequency analysis (|Laskaij|l993t ). we com- 
pute the main periods of the TTVs and the orbital elements. 
They are listed by decreasing amplitude of the trigonomet- 
ric term : the coefficient Ci represents the highest amplitude, 
and bold type ci indicates the precession rate of an angular 
variable in circulation. Let us note that only the period as- 
sociated with the trigonometric term of highest amplitude 
is indicated for 6i and to (the other period and/or a period 
combination are also present, as can be seen in Fig. [T|. 

Let us perform a similar study for a system of three 
planets in a Laplace resonance. We choose to analyze one 
of the Laplace-resonant co-planar syste ms formed by mi- 
:ratiori in a gaseous disc in the work of lLibert fc Tsiganij 
2011bl ). It consists of three planets with equal masses of 
1.5Mj {Mj — Jupiter's mass) and orbital periods of 20.1, 
40.5 and 81.2 days respectively (see Fig. ^. The TTVs of 
the inner planet are shown in Figure |3] and exhibits three 
periods now. Results of the frequency analysis are given in 
Table [1] (bottom). As in the two-planet resonant case, the 
long period of the TTVs signal corresponds to the preces- 
sion rate of the arguments of the pericenters. The other two 
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Figure 4. Comparison of the TTVs of the inner planet for the 
system given in Figure [2] (top panel) and for a modified system 
where the outer planet has been removed from the Laplace reso- 
nance (bottom panel). 



periods are related to 9i and 02, the libration angles of the 
2:1 mean-motion resonance between the two inner planets, 
and ^3 and O4,, the libration angles of the 2:1 mean- motion 
resonance between the two outer planets. These two periods 
are both easily detected in the time evolution of the Laplace 
resonant argument = Ai — 3A2 + 2A3 (see Figure [2] top 
right). 

The three periods in the TTVs signal are a signature of 
the multiple resonance, it is not due to the fact that the sys- 
tem consists of three planets. To ilustrate this, we compare 
in Figure S] the TTVs of the same three-body system as be- 
fore and another where the outer planet has been removed 
from the Laplace configuration. 



3 ON THE POSSIBLE DETECTION OF 
TERRESTRIAL PLANETS 

We now investigate if the three-periods signature described 
in the previous section survives when considering one or 
more Earth-like planets in the Laplace resonance. Three dif- 
ferent configurations are examined, depending on the num- 
ber of terrestrial planets in the system. 
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Figure 5. TTVs of the inner planet of the Laplace- resonant 
three-plan et Gliese 876 system . The parameters are issued from 
Table 3 of lRivera et al.l 1I2OIOI) . 



3.1 One terrestrial and two giant planets 

This case can be understood by studying Gliese 876 
(|Rivera et al.l (|201Ct l). It consists of two Jupiter- like plan- 
ets (0.7Mj and 2.3Mj) and a Uranus-mass body (14.6Me, 
Me = Earth's mass) in the Laplace resonance with orbital 
periods of 30.1, 60.1 and 124.3 days, respectively. The TTVs 
of the inner planet are given in Figure [5] Only two periods 
are detected in the curve, despite the Laplace configuration. 
This is confirmed with the frequency analysis. The same re- 
sult holds for a Jupiter- Jupiter-Earth system. 

Therefore, we conclude that the three-periods signa- 
ture cannot be detected for a Jupiter- Jupiter-Earth Laplace- 
resonant configuration. This is a consequence of the weaker 
gravitational effect of the smallest planet on its giant com- 
panions. Comparing Figures [S] and [T] we note that the TTVs 
curves look very similar. However, the presence of the third 
planet slightly modifies the period values and the amplitude 
of the signal. A deep analysis is needed to investigate if this 
difference can help to detect the third body. This is beyond 
the scope of this paper. 



3.2 Two terrestrial and one giant planets 

This case is particularly interesting. We analyze a giant 
planet of IMj in the Laplace resonance with two super- 
Earths of 10Mb, with orbital periods of 13.1, 26.2 and 52.5 
days, respectively. As in Section [51 the giant TTVs show 
three periods, see top left panel of Figure [6] The first one 
is the shortest (~ 300 days), arising from the 2:1 resonance 
between the two inner bodies. To identify the two other pe- 
riods, the time variations of the Laplace-resonant angle (j> 
and the pericenter argument 012 are given. We remark that 
(j) and the TTVs curve reach extrema at the same times. The 
period associated with the argument iU2 also appears in the 
TTVs, see top right panel of Figure [S] where we consider, 
for comparison, the reduced system of the two inner planets 
only (in the 2:1 resonance). 

As a result, the existence of two super-Earth compan- 
ions of a giant planet in a Laplace-resonant configuration 
can be inferred from the TTVs of the giant. The amplitude 
of the TTVs signal is high (~ 2.2 hours), and the time nec- 
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Figure 6. Left: TTVs and dynamical evolution of Laplace-resonant systems consisting of a giant planet and two terrestrial planets, on 
a time interval of ~ 100 years ; a giant of IMj with P = 13.1 days -|- two super-Earths of lOM^; on the top, a hot Jupiter of 0.5Mj with 
P = 5.1 days -|- two Earth-like planets of 2Me at the bottom. Right: Comparison with the systems formed by the two inner resonant 
planets only. 
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Figure 7. TTVs and dynamical evolution on ~ 400 years for a system of three Earth-like planets (1.5Mb, left), three super-Earth 
planets {15Me, middle) and three Jupiter-like planets (l.SAfj, right). Periods of the planets are 31.9, 63.8 and 127.5 days, respectively. 
The same initial orbital elements have been used for the three systems. 



essary to identify the three periods using frequency analysis 
(criterion ; twice the longer period of the signal) is approxi- 
mately 2000 transits, equivalent to ~ 70 years in that case. 

Another example is given at the bottom left panel 
of Figure |6l with a system of a hot Jupiter of 0.5Mj 
{P = 5.1 days) and two terrestrial planets of 2AfB (P — 10.2 
and 20.5 days, respectively). Since the terrestrial masses are 
smaller than in the previous case, ~ 3000 transits {^^ 40 
years) are needed to detect the three periods by frequency 
analysis. 



3.3 Three terrestrial planets 

Finally, let us consider a Laplace-resonant system of three 
terrestrial planets (1.5M_b) with periods of 31.9, 63.8 and 
127.5 days, respectively. The TTVs curve of the inner planet 
is shown in Figure [7] (left panel). As expected, the amplitude 
of the TTVs signal is extremely high (~ 3.5 days), but only 
two periods are detectable. The long period associated with 
the argument uj is longer than the time-span considered, 
and thus visible by increasing the mass values up to 0.5Afj 
(middle and right panels). Indeed the smaller the masses of 
the planets, the longer the precession rate of the arguments 
of the pericenters. 

Therefore, for a system of three Earth- like planets in a 
multiple resonance, the secular period cannot be detected 
in the TTVs in a time-span of 100 years. Only the two 
periods associated with the Laplace angle can be deduced, 
the longer being, from our example, ^ 1500 transits (~ 120 
years) for 1.5AfT planets. Nevertheless, the double-frequency 
TTVs signal is here the signature of a double 2:1 reso- 
nance between the three Earth-like bodies, with amplitudes 
of ~ 42 and 83 hours. This can be retri eved from two-bod y 
resonance calculations (e.g. Section 6 of lAgol et al.l (|2005h ). 
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Figure 8. Double period TTVs for different system configura- 
tions on ~ 40 years. Top panel: the three Jupiter-like planets 
given in Figure [7] Bottom panel: the Gliese 876 two-planet sys- 
tem shown in Figure [T] In each case the transiting giant planet is 
~ O.bMj with an orbital period of 30 days. 



© 2012 RAS, MNRAS 000, [T]-?? 



TTVs of Laplace-resonant systems 7 



4 DISCUSSION OF THE RESULTS 

4.1 Limitations due to observation time-span 

We have shown in the previous section that at least sev- 
eral decades are necessary to practically identify Laplace- 
resonant three-planet systems containing terrestrial bodies. 
On a given short observation time-span, only two periods 
can usually be detected in the TTVs of a Laplace-resonant 
system, thus leading to a curve very similar to the one of 
two planets in resonance. This is illustrated in Figure [8] 
the long-period resonant argument of a three- Jupiter sys- 
tem matches the pericenter precession of a two-body config- 
uration. This example suggests possible degeneracy between 
two- and three-planet resonant configurations on short ob- 
servation time-spans, and this question should be investi- 
gated in more details in the future. 

Therefore, long-term TTVs observations are necessary 
to constrain dynamically multiple-planet resonant systems, 
especially if no additional data such as radial velocity data 
are available. Let us note that, even if this work identifies the 
three periods in the TTVs signal as a signature of a resonant 
three-body configuration, the same three periods could cor- 
respond to different parametrizations of the system, as the 
degeneracy in the characterization of non-transi ting planets 
point ed out for a two-body resonance (see e.g. iBoue et al.l 
|2012| ). Furthermore, only coplanar systems have been in- 
vestigated here. The detection of inclined Laplace-resonant 
systems should also be adressed, as the possible formation 
of (highly) inclined two- and three-b ody resonant systems 
has recently been demonstrated (e.g . iThommes fc Lissaueij 
l2003l :[ Libert fc Tsiganis|[2009l . l2011bl ) . 



4.2 Kepler systems 

At the time of writing, three c onfirmed Kepler syst ems con- 
sist of three planets : Kep ler-9 (|Holman et al.l 2010l). Kepler- 
18 (ICochran et al.1 1201 ll ) and Kepler-30 (|Fabrvckv et all 
12013 ). The TTVs of Kepler-9b and Kepler-9c present the 
characteristic signature of a 2:1 mean-motion resonance. 
These two giant planets also have an inner non-resonant 
super-Earth companion. Kepler-18 exhibits a similar dy- 
namics, with two Neptune-mass planets in a 2:1 mean- 
motion resonance confirmed by TTVs and an inner non- 
resonant super-Earth planet. As argued in Section FS.ll such 
systems could host a terrestrial Laplace-resonant planet, but 
this would require further investigations. Kepler-30 is an in- 
teresting three-planet system with significant TTVs. While 
the Neptune-size Kepler-30b lies interior to a 2:1 mean- 
motion resonance with the Jupiter-size Kepler-30c, an addi- 
tional Jupiter-size planet orbits the star with a period that 
roughly corresponds to 5 times the one of Kepler-30b and 
2.5 times the one of Kepler-30c. The possible proximity to 
such a three-planet resonance could explain the difficult in- 
terpretation of the TTVs on the current short observation 
time. A longer time-span will help to identify the possible 
signature of three periods in the signal - if they exist. Let 
us note that most recent Kepler systems (e.g. Kepler-51, 
Kepler-53 and Kepler-60) sho w potential multibody reso- 
nances (see ISteffen et al]|2012l ). Furthermore a lot of KOIs 
are close to multi ple period ratio comm ensurabilities For in- 
stance, KOI-500 iLissauer et al.ir201ll ) is particularly inter- 
esting with possibly two three-body resonances. However, a 



study of the inner TTVs curve only is not sufficient for the 
dynamics of a five-body system. 



5 CONCLUSION 

We studied the possible detection of Laplace-resonant three- 
planet systems on a maximal time-span of ~ 100 years, as- 
suming that only the inner planet transits the star. We have 
shown that a three-body resonant configuration has a spe- 
cific TTVs shape: three periods are present, the long one 
associated with the precession of the arguments of the peri- 
centers, and the other two with the time evolution of the 
Laplace-resonant angle. Our criterion for detecting a signal 
is that the frequency is sampled completely, not that the 
amplitude of the Fourier component is large. 

The three periods are detectable for a giant planet per- 
turbed by two terrestrial bodies. A system of three Earth- 
like planets can be constrained from the inner planet TTVs, 
even if only two periods are visible. We discussed the num- 
ber of transits necessary to secure a detection. We pointed 
out that two and three-body resonant configurations could 
be mistaken on a short observation time, as the long-period 
resonant argument of a three- Jupiter system can be confused 
with the pericenter precession of a two-body configuration. 

Many Kepler extrasolar system s (and KOIs) reveal to 
be in a closely packed configuration (|Lissauer et al.l (|201ll )). 
and multiple-planet resonances can act as a phase-protection 
mechanism which insures long-term stability. Formation sce- 
narios (during the gas-disk phase) can also favor multiple 
resonant systems. We emphasize the importance of TTVs 
long-term observing programs for the detection of such sys- 
tems in the future and the contraints it will provide on the 
formation history of planetary systems. 



ACKNOWLEDGEMENT 

We thank the anonymous referee for his/her useful sugges- 
tions. The work of A-S L is supported by an F.R.S.-FNRS 
Postdoctoral Research Fellowship. A-S L thanks the people 
at Lille Observatory for their hospitality during her stay. 



REFERENCES 

Agol E., Steff'en J.H., Sari R., Clarkson W. et al., 2005, 

MNRAS, 359, 567 
Agol E., Steffen J.H., 2007, MNRAS, 374, 941 
Bean J.L., 2009, A&A, 506 

Borucki W. J., Koch D. G., Basri G., Batalha N. et al., 

2010, Science, 327, 977 

Borucki W. J., Koch D. G., Basri G., Batalha N. et al., 

2011, ApJ, 728, 117 

Bone G., Oshagh M., Montalto M., Santos N.C., 2012, MN- 
RAS Letters, 422, L57 

Cochran W.D., Fabrycky D.C., Torres G. et al., 2011, 
ApJS, 197, id. 7 

Collier C.A., Bouchy F., Hebrard G. et al., 2007, MNRAS, 
375, 951 

Csizmadia Sz., Renner S., Barge P., Agol E. et al., 2010, 
A&A, 510, id.A94 



© 2012 RAS, MNRAS 000, [l]-?? 



8 A.-S. Libert and S. Renner 



Fabrycky D.C., Ford E.B., StefTen J.H. et al., 2012, ApJ, 
750, id. 114 

Ford E. B., Howell M. J., 2007, ApJ, 664, L51 

Ford E.B., Rowe J.F., Fabrycky D.C., Carter J.A. et al., 

2011, ApJS, 197, Issuel, id. 2 

Ford E.B., Fabrycky D.C., Steffen J.H., Carter J.A. et al., 

2012, ApJ, 750, Issue 2, id. 113 

Holman M.J., Murray N.W., 2005, Science, 307, 1288 
Holman M.J., Fabrycky D.C., Ragozzine D., Ford E.B. et 

al., 2010, Science, 330, 51 
Kipping D.M., 2009, MNRAS, 392, 181 
Laskar J., 1993, Physica D, 67, 257 
Libert A.-S., Tsiganis K., 2009, A&A, 493, 677 
Libert A.-S., Tsiganis K., 2011a, MNRAS, 412, 2353 
Libert A.-S., Tsiganis K., 2011b, CeMDA, 111, 201 
Lissauer J. J., Ragozzine D., Fabrycky D.C. et al., 2011, 

ApJS, 197, id. 8 
Miller-Ricci E., Rowe J. F., Sasselov D., Matthews J. M. 

et al., 2008, ApJ, 682, 586 
Miralda-Escude J., 2002, ApJ, 564, 1019 
Morbidelli A., Tsiganis K., Crida A., Levison, H.F., Gomes 

R., 2007, AJ, 134, 1790 
Nesvorny D., Morbidelli A., 2008, ApJ, 688, 636 
Nesvorny D., Kipping D. M., Buchhave L. A., Bakos, G. . 

et al, 2012, Science, 336, 1133 
Reidemeister, M., Krivov, A.V., Schmidt, T.O.B. et al., 

2009, A&A, 503, 247 
Rivera E.J., Laughlin G., Butler R.P., Vogt S.S., Haghigh- 

ipour N., Meschiari S., 2010, ApJ, 719, 890 
Schneider J., Dedieu C, Le Sidaner P., Savalle R., Zolo- 

tukhin I. et al., A&A, 532, A79 
Steffen J.H., Agol E., 2005, MNRAS, 364, L96 
Steffen J.H., Fabrycky D.G., Ford, E.B., Carter J.A. et al., 

2012, MNRAS, 421, 2342 
Thommes E.W., Lissauer J.J., 2003, ApJ, 597, 566 



© 2012 RAS, MNRAS 000, [T]-?? 



